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Abstract 



A method for computing the stress-energy tensor for the quantized, massless, 
spin 7j field in a general static spherically symmetric spacetime is presented. 
The field can be in a zero temperature state or a non-zero temperature thermal 
state. An expression for the full renormalized stress-energy tensor is derived. 
It consists of a sum of two tensors both of which are conserved. One tensor 
is written in terms of the modes of the quantized field and has zero trace. 
In most cases it must be computed numerically. The other tensor does not 
explicitly depend on the modes and has a trace equal to the trace anomaly. It 
can be used as an analytic approximation for the stress-energy tensor and is 
equivalent to other approximations that have been made for the stress-energy 
tensor of the massless spin ^ field in static spherically symmetric spacetimes. 
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I. INTRODUCTION 



One of the most useful tools in the study of quantum effects in curved spacetime is the 
stress-energy tensor (T^,). The stress-energy tensor gives information about the effects the 
spacetime has on the quantized fields that propagate in it and also provides a source for 
that geometry via the semiclassical backreaction equations 

G„ v = 8tt(T^) . (1.1) 

Unlike alternative approaches, such as using Bogolubov transformations or particle detec- 
tors, the stress-energy tensor includes the effects of both particle production and vacuum 
polarization. It therefore does not depend on the distinction between these two effects which 
is often blurred in curved spacetime. 

In four-dimensions, the computation of (T^) for a quantized field in a given curved 
spacetime background is at best nontrivial and at worst almost impossible, given current 
technology. There are four types of approaches that have been used to compute the stress- 
energy tensor for a quantized field in a static spherically symmetric spacetime. One is to 
constrain the form of the stress-energy tensor by integrating the conservation equation and 
using the symmetry properties of the state that the field is in. If the field is conformally 
invariant the fact that its trace is known can also be used. This approach has been used in 
Schwarzschild spacetime |lj for fields in states that are static and spherically symmetric and 
for conformally invariant fields in such states in a global monopole spacetime |2||| . 

A second approach is to derive an analytic approximation for the stress-energy tensor. 
This approach is sometimes combined with the previous one. Analytic approximations have 
been derived for conformally invariant fields of spin 0, |, and 1 in Schwarzschild space- 
H~n]], in Reissner- Nordstrom and extreme Reissner- Nordstrom spacetimes 
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time 

and in a general static spacetime [|2U[| . An approximation has also been derived for mass- 
less scalar fields with arbitrary curvature coupling in a general static spherically symmetric 
spacetime pl |. The approximations in Schwarzschild and Reissner-Nordstrom spacetimes 



are mostly for the Hartle-Hawking-Israel f22| , |2"3|| state but in some cases are also for the Un- 
ruh |24|] and Boulware [ 25| states. The more general approximations are for the vacuum state 
and a thermal state at an arbitrary temperature. For massive fields the DeWitt-Schwinger 
expansion |2(J can often be used to obtain a good approximation for the stress-energy ten- 
sor. This approximation has been derived for massive scalar fields in Schwarzschild space- 
time [p^ -|30|, in a general static spherically symmetric spacetime [ 2l] , |l? |, and in a general 
spacetime [ 31 ] . It has also been derived for massive spin | and 1 fields in Reissner-Nordstrom 
spacetimes ||32| . The DeWitt-Schwinger approximation is independent of the state of the 
field. 

A third approach is to numerically compute the stress-energy tensor. This has been 
done for conformally invariant spin and spin 1 fields in Schwarzschild spacetime on the 
event horizon |33fl and outside of the event horizon |R^|5l,|35l|2Tl,|36| . Numerical calculations of 



the stress-energy tensor for both massive and massless scalar fields with arbitrary curvature 
couplings in Schwarzschild and Reissner-Nordstrom spacetimes outside the event horizon 
have also been done [21,0]- Finally a numerical computation of the stress-energy tensor for 
the massless spin | field has been done in a global monopole spacetime [38 . 

The fourth approach is the most accurate and the most difficult. It is to analytically 
compute the full renormalized stress-energy tensor. This has been done for the special case 
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of conformally invariant fields in conformally flat spacetimes [|39| , f40| and for scalar fields with 
arbitrary mass and curvature coupling in de Sitter space 1 41- 45)1 . The stress-energy tensor 
for the massless spin 1 field has also been computed analytically on the event horizon of a 



Schwarzschild black hole [}4q] . To our knowledge no other exact analytic calculations of the 
stress-energy tensor in static spherically symmetric spacetimes have been published. 

Clearly only a little has been done regarding the computation of the full renormalized 
stress-energy tensor for the massless spin | field. In particular the only calculation of the full 
renormalized stress-energy tensor in a black hole spacetime is on the event horizon of the ex- 
treme Reissner- Nordstrom black hole [ |T7[] . This was possible because the geometry near the 
horizon asymptotically approaches the conformally flat Bertotti-Robinson geometry f4"T|-|4"9 



In this paper we present a practical method that can be used to numerically compute the 
stress-energy tensor for the massless spin | field in a general static spherically symmetric 
spacetime. The method makes use of the Euclidean Green's function for this field. Our 
approach closely follows that of Anderson, Hiscock, and Samuel [21| which in turn was built 
on that of Howard and Candelas []5]J| . As in those papers, the Euclidean Green's function 
is constructed using a method similar to that of Candelas P^| , and point splitting is used 
to renormalize the stress-energy tensor. The point splitting counterterms for spin 0, |, and 
1 fields have been computed by Christensen [EnjEjlj]. 

Our main result is the derivation of a renormalized expression for the stress-energy tensor 
of the massless spin | field in a general static spherically symmetric spacetime. The field can 
be either in a vacuum state or in a thermal state at an arbitrary temperature. One advantage 
of our expression is that one does not need to work with a vierbein or any Dirac matrices 
to use it. These were involved in the derivation of the expression but no longer are a part 
of it. In our derivation the stress-energy tensor is divided into two parts. One part depends 
on sums and integrals over radial mode functions. This part usually needs to be computed 
numerically. The other part consists of an analytic tensor with a trace equal to the trace 
anomaly for the spin | field. Since both parts are separately conserved, the analytic part 
can be used as an analytic approximation for the stress-energy tensor. It is equivalent to 
the approximation of Brown, Ottewill, and Page @ in Schwarzschild spacetime, and to that 
of Frolov and Zel'nikov |2(J in a general static spherically symmetric spacetime for certain 
values of the three arbitrary coefficients in their expressions. 

In Section [H], we review the formalism used to describe spinors in curved space and 
its application to static spherically symmetric spacetimes. In Section |TJ we construct an 
expression for the Euclidean Green's function for the massless spin | field in a general static 
spherically symmetric spacetime. The WKB approximation for the radial mode functions 
is introduced in Section [TV] and the point splitting renormalization procedure is reviewed 
in Section |V[ In Section [Vl| an expression for the unrenormalized stress-energy tensor is 
derived. The renormalized stress-energy tensor is derived in Section [VI H , and a part that 
can be used as an analytic approximation is displayed. Our results are summarized and 
discussed in Section |VIII| . Various technical results are derived in the Appendices. 

Our conventions are those of Misner, Thorne, and Wheeler and we use units such 
that h = c = G = k B = l. Greek indices are used to refer to general spacetime coordinates 
and Latin indices refer to locally inertial coordinates. 
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II. REVIEW OF SPINORS IN CURVED SPACETIME 



The study of spin | fields in curved spacetime appears to have begun with Fock and 
Ivanenko |53| who derived the Dirac equation in a general spacetime. Brill and Wheeler [54 



adapted the formalism of spinors in curved space to static spherically symmetric spacetimes 
and derived the specific form of the Dirac equation in these spacetimes. Green's functions 



for the spin | field have been derived by Boulware [55] for Schwarzschild spacetime and by 
Bezerra de Mello, Bezerra, and Khusnutdinov |38j in a global monopole spacetime. 

In this section we review the formalism most commonly used to treat spin | fields in 
a general curved spacetime. We then derive an explicit form for the Dirac equation in 
Euclidean space and discuss its solutions. 



A. General Results 



The spin of a field ip is defined by its behavior under a Lorentz transformation. An 
infinitesimal Lorentz transformation in flat space can be written in the form 



A 



a _ ro 



b + U) b 



(2.1) 



with u a b an arbitrary infinitesimal antisymmetric tensor. The transformation law for a field 
ip of arbitrary spin is |)B 



V -> D(A)if) = exp (\u ab a ab ) $ 



(2.2) 



where D denotes a representation of the Lorentz group and the matrices a ab are the gener- 
ators of this representation. For an infinitesimal Lorentz transformation, 



D(A) = 1 + \uj ab & 



ab 



The matrices a ab satisfy the commutation relations 



_a6 „cd 

a , a 



^cb^ad _ ^ca^bd _|_ ^db^ca _ ^da^cb 



Here rj ab is the usual Minkowski metric in Cartesian coordinates. For spin \ fields 



a 



ab 



where the 7° are 4x4 matrices satisfying the anticommutation relations^ 

{Y,i b } = -2 V ab i. 

For this work, we use the Dirac representation [|57j[: 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



1 Had we chosen {j a ,"f b } = 2r] ab l, we would have chosen a a 
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y=t; _»), y = (_° ( ;-j -1,2,3. (27) 

Here the <7j are the usual Pauli matrices 



In curved spacetime, it is useful to use the vierbein or tetrad formalism |p6[ . In this 
formalism, one sets up a locally inertial frame at every point. The coordinate transforma- 
tion between the general spacetime coordinates and the coordinates for the frame that is 
locally inertial at a given point x is given by the vierbein functions e a ^(x). As a result the 
relationship between the locally inertial metric 7] ab and the spacetime metric g^ v at x is 

Vab = g^u{x)e a ll (x)e b v (x) . (2.9) 

Note that the locally inertial coordinate indices are raised and lowered with 7] ab and the 
general coordinate indices with g^ v . It is easy to show that the vierbein functions also 
satisfy 

e/e\ = 5». (2.10) 
Using the vierbein formalism, one is able to maintain Lorentz covariance and general 



coordinate covariance through the definition of an appropriate covariant derivative [p6 
One accomplishes this by requiring that the spinor covariant derivatives behave in the same 
manner in curved space as in flat space. In flat space the gradient of a spinor is both a spinor 
and a vector. This implies that under a Lorentz transformation d^ij) — > A a ^D (A) dpip. In 
curved space, the covariant derivative is defined in such a way that this relationship is 
preserved |56| , |58| . Thus 



V a ij(x) -> A a b (a;) D (A (x)) V b ^{x) (2.11) 

where V a = e a M V M and A is a function of x, since the locally inertial frame is different, in 
general, at different points in the spacetime. Equation Q2.ll ) is satisfied if we define 



V M = <9 M + IV (2.12) 

with 

T M = \a ab e a v e bv . tlJL . (2.13) 

Here and throughout the semicolon denotes the covariant derivative of a vector or tensor. 
It acts on the spacetime indices (denoted by Greek letters). 

With the vierbein formalism, the Dirac equation for a massless spin | field in an arbitrary 



spacetime can be written as J58 



tYV^ = , (2.14) 
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where 7 M = 7 a e a M . 

To calculate (T^) in a static spacetime, it is advantageous to work in Euclidean space 
instead of Lorentzian space. To compute covariant derivatives in Euclidean space, one may 
repeat the above steps, but this time preserving covariance under a Euclidean rotation 
instead of a Lorentz transformation. An infinitesimal rotation TZ ab in Euclidean space can 
be written as 



K, 



ab 



Sab + ^ b (2.15) 

where u® b is an arbitrary antisymmetric tensor. Under a rotation TZ, a field ip of arbitrary 
spin has the transformation 



ab u E 



■0 — > D('R.)c = exp ( a 
For an infinitesimal rotation this becomes 



D(7e) = 1 + b»*o? 



(2.16) 



(2-17) 



The crff must be generators of a representation of the Euclidean group (in four dimensions). 
They satisfy the commutation relations 



E i°E 



£cb ad 

6 a E 



rca bd i zdb ^ 

o a F + o a 



db ^ca 
E 



rda cb 

6 a E 



For spin ~ fields 



with 



ab 



1e,Te 



(2.18) 



(2.19) 



{7^,7s} = 25 ab . (2.20) 
It is easy to show that this anticommutation relation is satisfied if we make the identification 



ie = -n 



(2-21) 



Note that the relationship between the vierbein and the metric in Eq. ( |2.9| ) allows us to 
use the same vierbein for both the Euclidean and Lorentzian sectors of a given static space- 
time. The Euclidean covariant derivative can then be defined analogously to the Lorentzian 
one, by requiring that the covariant derivative of a spinor in curved space behave in the 
same manner as in flat space: 



Vj^ (z) 
Eq. ( ggg ) is satisfied if 



nj (x)D(-Jl (aO)VfV (z) 



E 



with 



r 



2° E e a 



e-bv 



■41 



(2.22) 



(2.23) 



(2.24) 



Then the Dirac equation for a massless spin 1 field in an arbitrary Euclidean spacetime 



is 



- W = • 



(2.25) 
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B. Static Spherically Symmetric Spacetimes 



For a general static spherically symmetric spacetime the metric can be written in the 
form 



ds 2 = -f(r) dt 2 + h(r) dr 2 + r 2 d6 2 + r 2 sin 2 (#) 



^2.26) 



The transformation to Euclidean space is made by setting t = —it in Eq. (|2.26|) . One 
vierbein that is appropriate for this metric is 

/ f- 1/2 

/i _1 / 2 sin#cos0 r _1 cos#cos0 — r _1 csc sin i 

h^ 1 / 2 sin 9 sin r _1 cos#sin0 r _1 esc 9 cose/ 

V h~ 1 / 2 cos9 -r- 1 sin9 



^2.27) 



/ 



The rows of this matrix are specified by the index a and the columns by the index /i. One 
may derive explicit forms for the spin connection matrices by substituting Eqs. ( |2.5| ) and 
( [2 .27| ) into Eq. (|2.13|) . The spin connections for Euclidean space can be found by substituting 
Eqs. (|2~T9| ), ( ggg ), and (|2~27j ) into Eq. (ggg ). In both cases the results are 



r t 



IT 

o 



f7°7 



with 



4/1/2/^1/2 

|(1 - h- l/2 ) ( 7 1 7 3 cos0 + 7 2 7 3 sin0) 

|(/t -1 ^ 2 — 1) ^7 X 7 2 sin 9 + 7 1 7 3 cos 9 sin — 7 2 7 3 cos 9 cos0^ sin 6 1 
7 = sin # cos 7 1 + sin 6 1 sin 7 2 + cos # 7 s . 



^2.28) 



(2.29) 



Substituting Eqs. ( |2.28| ) into Eq. ( |2.25| ) results in the following form for the Euclidean space 
Dirac equation: 



where in the Dirac representation, 



a 
a 



(2.30) 



f2.31} 



and L is the standard angular momentum operator with components 

Li = i sin (9,9 + i cot # cos cL 



£2 
£3 



— i cos 4>dg + i cot # sin 



(2.32) 



Using separation of variables, one finds that, as in the Lorentzian case [ pl| , p9] , |50| , there 
are two types of solutions to the Euclidean Dirac equation. When this equation is written 
in the form (|2.30|) , they are 
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F wJ (r)tf£ + (0,. 



r / 1/4 \G Utj {r)V£_ 
e 1 ^ {G Uij (r)V£_ 



r fl/4 



[2.33) 



The angular functions ty™± are eigenf unctions of the operator k = a ■ L + 1, and satisfy 



\&™ ± = ±(j + \) tt£ ± and crtf^ 



where 



^2.34) 



cr = sin 9 cos oi + sin 6* sin rr 2 + cos # cr 3 . 



^2.35) 



Their explicit form is 



J 2 

1 



3-m y m 1 2 



2.7 



/ l j+l-m y m -\ \ 



and 



2j+2 



3 + o 



V 



i±l±my"' T : 
2i+2 • , 1 



(2.36) 



with j 



I 3 

2' 2' 



and m 



-j, — j + j. The Yj ^jj 2 are ^ ne standard spherical harmonics. 
Substituting either of t/>i(o;) or "02 (#) in Eq. (|2.33|) into Eq. (|2.3CI|) yields the following 
coupled set of mode equations for F U j and G^f 



drG^.j H Cr. 



„ 1 „ „ j+i 



-d F — 



fl/2 ^ h 1 / 2 / 

These may be combined to form the second order equations 



(2.37a) 
(2.37b) 



h l/2°r 
fl/2 



fl/2 

t f) f . 



h 1 /' 



f 



1/2 



^^drGujj 



fl/2 
r ">3 fol/2 Ur 

fl/2 

+ G w j , 1 /r> d r 



a/ 



1/2 



A 2 / 



/lV2 



- H) . 
/ 



2 u>3 



u 2 G 







(2.38a) 



(2.38b) 



Since F U j and G^j obey a coupled set of linear first order equations, there will be two linearly 
independent solutions. Boundary conditions for the radial modes depend on the physical 
situation and are usually specified at both large and small values of r. For example, in a 
black hole spacetime the modes must be finite at the event horizon and at r = oo. We 
denote the solutions that satisfy the boundary condition at small r by the superscript p. 
These are F% j and G^j- Those that satisfy the boundary condition at large r are denoted 
by the superscript q. These are F* ,- and Cr^,,-. The normalization of these solutions is, for 
the moment, left unspecified. 
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III. THE EUCLIDEAN GREEN'S FUNCTION 



The Euclidean Green's function for a massless spin ~ field is required to be a solution of 
the equation 



-j E V°S E (x,x') 



8 A [x — x') 



(3-1) 



It must also satisfy the appropriate boundary conditions. 



Using a method similar to that which Candelas |33[ used to construct the Euclidean 
Green's function for a scalar field, we find that the Euclidean Green's function for the spin 
i field in a zero temperature vacuum state in a static spherically symmetric spacetime can 
be written in the form 



S e (t, x; t', x 1 



(r—r 



F^(r)^ + (9, 



duo ooe^-^ ~ * 
2 7 rrr7(r) 1 /4/( r ')i/4 2^ 



'Wr m t //V J.A 



+ 0{r'-r) 



(3.2) 



The notation ® denotes a direct product. The notation p <-> g means that the superscripts 
p and g on the radial mode functions F u> j and G u j should be switched, but everything else, 
including the arguments of the mode functions, should be left the same as in the previous 
expression. Thus the explicit form of the second term in Eq. (|3.2j ) is 



r — r 



■f*,(r)*£ + (M) 



® (G^(0*,, 



For a thermal state at temperature T the integral in Eq. 
the following rule: 



is changed to a sum using 



/oo °° 
duoV{uo) -^k ^ [V{tik) + V(-tik)] (3.3) 

■°° n=l/2 

where k = 2nT and the sum over n runs over the half integers |, |, . . . . This ensures that 
the spin | Euclidean Green's function is antiperiodic in Euclidean time r as it must be for 
an anticommuting field | J58|| . 

An examination of Eq. ( |3.2| ) shows that it satisfies the same boundary conditions that 
the radial modes do at large and small values of r. To verify that Eq. (|3.2|) can satisfy Eq. 



9 



(|3.1|), first note that because Eq. ( |2.33|) is a solution of the homogenous Dirac equation, most 
of the terms in ^^V^Se will vanish. The only non-vanishing terms will be those where the 
derivative acts on the radial step-functions, yielding factors of S(r — r'). Combining terms, 
and writing out the matrices explicitly, many of the terms will vanish. If the condition 



GWr)f*j(r) - F^Gl^r) ] = 1 (3.4) 



is imposed, then the integral over uj can be computed to yield a factor of 5(t—t'). However, 
Eq. fl3.4p is just a Wronskian condition for the solutions to Eqs. ( |2.37a| ) and ( |2.37b| ), so the 
normalization can be chosen so that it is satisfied for all values of r. All that is left is the 
angular mode sums, which can be computed using the identity 

£ £ me, wne; 00 = »«> - - m , (3 . 5 ) 

e=om=-e sm6 



and the proof that Eq. ( |3.2| ) satisfies Eq. (|3.1|) is complete. 

In spite of the fact that Eq. ( |3.2|) is a solution to (|3.1|) with the correct boundary con- 
ditions, it might appear that our ansatz is not sensible. After all, in Lorentzian space, we 
expect expressions like ip(x)ip(x'), and not ip(x)if> '(x 1 ) to appear in the Green's function. 
For example, one knows that ipif> is invariant under Lorentz transformations. Since the trace 
of our Se when the points are brought together contains if) 'if) instead of ipip = ifj^'y ip, at 
first glance it may appear unreasonable. 

The reason for the absence of the 7 is that we are working in Euclidean, not Lorentzian 
space. A finite rotation in Euclidean space of a field if) will transform according to Eq. ( [2.161) . 



The a% s in Eq. (|2.16|) are all anti-Hermitian, as can be seen easily from their definition in 
Eq. ( gig) . It follows that 



- ^ exp (iu&rjt) = ^ exp (-§u;>| b ) , (3.6) 

so that under a local Euclidean rotation, if)<if) — > ip^ip. 

In contrast, in Lorentzian space, the cx afe 's are not all anti-Hermitian. Specifically, 



a 



°*t = a-oi while a ut = _ a ij . (3.7) 



Here i and j are space indices. It is easy to demonstrate that cr a ^7° = — 7°cr ab , so that 
under a Lorentz transformation, 

^t 7 o ^ ^ exp (l^a^) 7 ° = exp (-|o; o6 cr a6 ) (3.8) 

and thus ip^°ip ip^^ip. Therefore the overall form of our Se as the direct product 
if>(x) <8) if)^{x') is sensible. In addition, we show in Appendix [A] that our Green's function is 
equal to the standard one in flat space. 

It is possible to compute the sum over m in Eq. ( |3.2|) by using the addition theorem for 
spherical harmonics fBTfl: 



4tt 1 

^(cos 7 ) = e Yne^yno'A') , (3.9) 
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with 



The result is 

S e (t, x; t',x') 



cos 7 = cos 9 cos 9' + sin 9 sin 9' cos(</>— <//) 
d\xuj 



(3.10) 



x 



X 



rr'(f(r)f(r')y/ 4 
J2\e(r-r'){-sm[u;(T-T')} 7° A 9 J p e + i cos[lo(t-t')] 7^} 



£=0 



Pticos'y) 



+ 9(r'—r)^p <-> gj 

where / d\i is defined such that for an arbitrary function V{ui) 

[ da V(u) = At r duoViuo) if T = 
J An 2 Jo 

if T > 



(3.11) 



4^ D 

n=l/2 



(3.12) 



and 



2 ' 1 2 



?+!)! + £ -L 



A:;./ = -g „ , 1 (r)R 



(£+1)1 + 2- L +F«_ 1 (r)Gl e _ 1 _(r') 



1 — E-L 



(3.13) 



Note that the subscripts £ + | on the radial mode functions denote the value of j in the 
mode equations ( |2.37a| ) and ( |2.37b| ). The operators L are defined in Eq. ( [2.32j ), and always 
act on the unprimed variables 9 and 0. 



IV. WKB APPROXIMATION 



The method presented in Sections |\T] and |VII| of computing the stress-energy tensor for 
the massless spin | field makes use of a WKB approximation for the solutions to the mode 
equations (|2.37a|) and (|2.37b|) . The WKB approximation can be obtained by writing 



2W(r) 
1 

2W{r) 



: exp 



: exp 



W(r) (h/f) 1/2 dr 
W{r) (h/f) 1/2 dr 



Substitution of Eqs. (pLTaD and (pTTbT) into Eq. (|2.37a|) gives 



F q .(r) 



to 

CO 



fi/2 W'(r) , f 1 ^ 



W(r) 



h 1 / 2 2W{r) 
f 1 ' 2 W'{r) 
W 2 2W{r) ~ 



(3+ 



1/2 < 



(4.1a) 
(4.1b) 



(4.2a) 
(4.2b) 
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Note that Eqs. (|4.1a| -b) and (|4.2a| -b) are normalized in this form so that the Wronskian 
condition (|3.4j) is satisfied. Then substituting Eqs. ( |4.1a| ) and (|4.2a| ) into Eq. ( |2.38b| ) or 
substituting Eqs. ( |4.1bjj and ( [4.2b| ) into Eq. ( |2.38b| ) gives the following equation for W(r): 



W 2 (r) = n 2 (r) + 



fW(r)" 3fW(r)' 2 fh'Wjr)' f'W(r)' 
2hW(r) ~ AhW{r) 2 ~ Ah 2 W(r) + AhW{r) 



with 



Q 2 (r) 



2 / i l/2 r h l/2 r 2 



Eq. ( |4.3| ) can be solved iteratively. The zeroth order solution is W(r) 
order solution is 



W{r) = Q + 



4 



f 1 

1 n" + 



!!__!]£ 

2tt 2 [ h h? 
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2htt 3 



(4.3) 
(4.4) 

O(r). The second 
(4.5) 



and so forth. 



V. POINT SPLITTING 

As is the case for most quantities in quantum field theory, the stress-energy tensor for 
quantized fields is divergent, and must be renormalized. The method of renormalization 
which we choose is that of point splitting. In this method PB| one expands the expression 
for the stress-energy tensor in powers of the geodetic interval a(x, x') and its derivatives 
a' ,u (x,x'). The quantity a(x, x') is defined as one-half of the square of the proper distance 
along the shortest geodesic that joins the points x and x'. It has the property 

a(x,x') = l(a^) (5.1) 

with 

<r" = a'^ . (5.2) 

One then subtracts the results from the unrenormalized stress-energy tensor and allows the 
points to come together. Christensen [ pO| , pT| has calculated the renormalization counterterms 
for the stress-energy tensor for free quantized fields of spin 0, |, and 1. To actually implement 
the point splitting renormalization scheme, one usually chooses a set of coordinates, and 
chooses the way in which the points are split. The geodetic interval a and its derivatives 
are then expanded in powers of x M — 

As is shown in Section |V]|, the unrenormalized stress-energy tensor may be written in 
terms of covariant derivatives of the Euclidean Green's function. It is made up of terms 
which are bi-spinors and sometimes also bi-vectors. A bi-spinor is an object that transforms 
like a spinor at x and a spinor at x' under Lorentz transformations (or Euclidean rotations) 
and a bi- vector is an object that transforms like a vector at x and a vector at x'. 

In general, a spinor at x' is parallel transported to a spinor at x through the use of the 
bi-spinor of parallel transport I(x, x'). A vector at x' is parallel transported to a vector at x 
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through use of the bi- vector of parallel transport g^y ne may show that an object which 
is both a vector and a spinor is parallel transported by y I(x,x f ). For example, 

(7^(a0)|,= gVix^WW) (5.3) 

is a spinor and vector at x, and a scalar at x' . We show in Appendix [B] that to paral- 
lel transport tp(x') to ip(x), one must multiply on the right hand side by I(x',x), where 
I(x,x')I(x',x) = 1. Then 

ty(x)-f) y = g\4{x')l X ' I{x\ x) (5.4) 

is again a spinor and vector at x and a scalar at x' . 

We shall split the points along the time direction so that 

e = t-t\ r = r', 9 = 9', <\> = <\>' . (5.5) 

The quantities a^ 1 can then be expanded in powers of e. The result for a general static 
spherically symmetric spacetime is |2l| 



t_ / /2 3 If f 4 3 f' 3 h' 3f' 2 f"\ 5 7 
a_e+ 24//i e 120 \8f 2 h 2 + 16 //i 3 8//i 2 / e+ ^' 

a =-4h € "2i(-W + 4^J C + ° (e) 

a e = a ^ = . (5.6) 
The bi- vectors of parallel transport can also be expanded in powers of e with the result |21j 

g u ' = --- J^t 2 + ( — f — f -^- + f3k ' ) e 4 + 0(e 6 ) 

y / 8/ 2 /i ^384/ 3 /i 2 9Qf 2 h 2 192/ 2 /i 3 / 

^ = -n^' = -JL e - ( f ' 3 + f ' f " - f ' 2h ' ) e 3 + 0(e 5 ) 
y y 2fh \96f 2 h 2 A8fh 2 96fh 3 J 1 ; 

„< _ i r 2 ( /* rr \ 4 6 

5 h %fh? e 1 384/ 2 /i 3 96//i 3 192//i 4 r {e ' 



1 

1 



y ~ -^Tn ■ ( 5 - 7 ) 

r z sm # 

All other components vanish. In Appendix [B], we show that I{x',x) has the following 
expansion in terms of e: 

I{x',x) =/ 1 l + 7 7 / 2 , (5.8) 

with 

1 16//i 2 [ 256 / 2 /i 2 32//i 3 16//i 2 J 24 16 j ' 

• r ^^75 + ^^)^ + 0^) ■ 



4/1/2^1/2 I 6 4 /S/2^3/2 32/1/2^5/2 igJl/2^3/2 I g 
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VI. AN UNRENORMALIZED EXPRESSION FOR (T^) 

The stress-energy tensor for the classical spin | field is 



T^(x) = -- 



(6.1) 



where parentheses around tensor indices indicate symmetrization. For example = 
| (A^ + A UfJi ). In Appendix [E| we show that, for the massless spin | field in a general static 
spherically symmetric spacetime, 

(T^(x)) unrcn = - J hm Im {Tr [ 7(/i (v„) + 5|] - <?,) A 'V V [S E + S C E }) I(x', x)] } , (6.2) 

where Se is understood to stand for Se(x,x'). In Appendix [D] the charge conjugation 
operation is defined for a matrix operator U as 

U c = C(TJ) T C j , (6.3) 

where T means transpose and C is the charge conjugation matrix, given by C = i^j in 
the Dirac representation. 

One can write the components of (T^ v ) unren explicitly in terms of mode functions by 
substituting Eq. ( |3.11| ) into Eq. fl6.2|). In Appendix [D|, it is shown that when the points 
are split in the time direction it is possible to write the results in terms of scalar functions 
multiplied by various combinations of Dirac matrices.^ This makes it possible to compute 
the trace in Eq. ( |6.2| ) without explicitly solving the mode equations or computing the mode 
sums. The nonzero components of the unrenormalized stress-energy tensor are 



2Re< I 



-uj 2 cos[cu(r-r')]/i (g tt +g tt )A 1 -iuj sm[u(T-r')]I ig tr —A 3 



^ 2 
f 1 



-UJ 2 cos[u;(r-r')]I 2 g tr "-^A 1 - iusHu(T-T'))I 2 (g*-g«) jjjj^A 

L 
4/ 



+ucoslu}{T-T')]I 2 g tr "- — A 2 + iujsm[uj(T-T')]I ig tr ' 1 



+ £\ A, 



(Tr r )r 



2Re< 



+ iuj sin Lmt— r 



')]h (<7 rr -< 



A 3 - uj 2 cos[u (t -r' )]I 2 g rt ' A ± 



2 Since our 7^ has the opposite sign as that in Ref. [58], this introduces an apparent minus sign 



between our Eq. (fy) and Eq. (3.191) of Ref. |8|. 



3 It is actually possible to do this for an arbitrary separation of the points but the result is more 
complicated. 
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-u 2 C os[uj{r-r')]h (g rr +g rr ') —^A l + uco8[u(r-r')]I 1 (g rr +g rr ') —A 



+ iu sin [tu(r-r')]/i5f rt ' 

(Tq )unren {^4> )unren 2Re 

+ iuj sin[a;(r— r')]I 2 



r 



4/W/2 



Ai 



r 3 yi/2 
1 



— u; cos[c<j(r— r')]/ 1 y4 2 — ia> sin[uj (t—t')]I 2 A 3 



i Mi 



(6.4) 



where 



(£+l)F g l (r)F v ^-IGP 1 (r)G p t (r) + 



l- 2 



fl/2 



A 2 

A 3 



E 

£=0 

oo 

E 



£(£+1) ^r)** !(r) + F 9 l (r)G p x (r) 



?+l) r 2 a; 
+ ~27 



£(£+1) F 9 x (r)F p jM+G 9 1 (r)G p x (r) 



r 



J >*+2 



,e-i 



+ 



2J1/2 2/ 1 /2/ i l/2 



r 2 /' 



4p/2 h l/2 



(6.5) 



As in the scalar field case p3|,p|,P, |62| , |2Tl| , there is a superficial divergence for large £, which 
arises because the points are split along the time direction. As discussed in f5j, it cannot 
be physical since one knows that the Green's function is finite with the points split. Since 
5(t — r') = if the points are split, one can add multiples of it and its derivatives to 
Eqs. (6.4). In particular one can add terms that will cancel the large I divergences in these 
equations. In Eqs. fl6.5|) , the terms that do not contain mode functions are these additional 
terms. 



VII. RENORMALIZED EXPRESSION FOR THE STRESS-ENERGY TENSOR 

As mentioned previously, the renormalized stress-energy tensor is obtained by subtracting 
off the point splitting counterterms and letting the points come together. Thus 

(Tfiiz/ren = lj m ((7]^)unren — (T^ds) • (7.1) 
x'— >x 

The DeWitt-Schwinger point splitting counterterms, (T^ U ) DS have been computed by Chris- 
tensen [|51||. In his derivation he eliminates odd powers of o from the point splitting coun- 



terterms. This results in the cancelation of terms with odd powers of a p in the numerator 
since there are only even powers in the denominators of all terms. However, the elimination 
of these terms results in the elimination of finite contributions to the stress-energy tensor, 
which must be included if it is to be conserved. 
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The linearly divergent terms can be derived through the use of equations (2.19) and 
(3.11) in Ref. g|. We fincf] 



(iu/BS,lm 



48vr 2 



(— R a (ti-,u) + \g^R' a + R 



ti Ft- 



P;i 



- [AR ^ - 2R ;Qjt + 2R x{fl ) ]2 -2[R 7 ^ + R 1 ^ ) 



(7.2) 



If the mode equation could be solved exactly for a general static spherically symmetric 
spacetime, and the mode sums and integrals in Eqs. (|6.4|) and (|5.5|) could be computed 



analytically, it would be straightforward to compute (T M y) ren . However, for most static 
spherically symmetric spacetimes, this is not the case. As a result, the limit in Eq. ( |7. 1| ) is 
difficult to evaluate since, without a solution to the mode equations, one doesn't know how 



to expand (T^) um . en in powers of x^ 



x 



in 



To get around this problem we use the WKB 



approximation to isolate the divergences in (T M ^) unrcn . 
in detail in Ref. [Ell]. We shall briefly review it here. 



The procedure to do this is discussed 



First one substitutes Eqs. (|4.1a|) , ( |4.1b|) , ( f4.2a|) , and ( f4.2b| ) into Eqs.( |6.4|) and computes 
the fourth order WKB approximation for the variable W . This can be done by iterating Eq. 
731) twice. Then the sums over I are computed using the Plana sum formula 



00 f°° t°° dt 

J2g(j) = ±g{k)+ g(r)dr + t -^ r — [g(k + it) - g(k - it)) 

j—h J k JO C 1 



(7.3) 



and the results are expanded in inverse powers of to. Only terms are retained that contribute 
terms of order u~ l or less to the stress-energy tensor, since these are the ones that lead to 
ultraviolet divergences when the points come together. If the field is in a zero temperature 
state, a lower limit cutoff A is used in the integral over uo for terms that are proportional to 
u^ 1 . In this way the quantity we call (T AiJy ) WK Bdiv is computed. We can add and subtract 
this quantity from the right hand side of Eq. ( fr"7i~l) with the result that 



(-^u/)ren ( (-^"/lm// unren - (^)\VKBdiv) + ((Tf.u) 
numeric (^(lu ) analytic • 



WKBdiv 



(T^ v )ds) 



(7.4) 



Because all of the divergences in (T pu ) unren and (T m ^)ds are contained in (T Mi ,)wKBdiv, both 
(Tpv) numeric and (T Ati/ ) ana i ytic are finite. The tensor (T^) numeric must usually be computed nu- 
merically, while (T^) analytic uiay be computed analytically for an arbitrary static spherically 
symmetric spacetime. 

To compute (T^J) analytic one can use the formulas in Appendix [FJto evaluate the integrals 
over to in (T^,) W KBdiv Then the analytic continuation discussed in Appendix |B| is used 
to analytically continue to Lorentzian space. The details of this process are discussed in 



Since we are working with a complex spin ^ field rather than the real field used in Ref. [51], we 
effectively have twice as many fields. Thus we have multiplied the counterterms in Ref. |51] by 2. 
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Appendix |E[ The resulting terms are proportional to various powers of (t — if) or to ln(t — t'). 
The point splitting counterterms are then also expanded in powers of (t — t'). After the 
subtraction the limit f — > Ms computed. As mentioned above, for the zero temperature case, 
there is an infrared divergence in the a; -1 term which can be handled with the introduction 
of an infrared cutoff, A. To compute numeric one can simply set r' = r in the mode 
integrals since this quantity is finite in the limit that the points come together. 
Our results are: 



\T[/,ii/ analytic ('//jf)o (-^~/w)log 



(7.5) 



with 



1 



240vr 2 
1 



7k 4 
4f 



+ 



+ 



+ 



5k 2 
6r 2 f 

A3f' 2 



5k 2 



5k 2 /' hn 2 h' 

I — H 

Qr 2 fh 3rf 2 h 6rfh 2 

5f 2 19f 3 



~ 12f 2 h 2 
77/' 4 



+ 



24/ 3 /i 6/ 2 /i 
2/i' 15f'h' 



h 2 2Ar 2 f 2 h 2 8r 2 f 2 h \2rf 3 h 2 192ph 2 



bf'ti lhf' 2 ti f A ti 



f/3l 



12r 2 fh 2 
5/" 



+ 



7/i 



/2 



/2 



+ 



6r 2 fh 2 
h" 

+ 



+ 



8rf 2 h 3 2Af 3 h 3 4r 2 /i 4 
5/" 5/'/" fV" 



+ 



2r//i 4 
2h'f" 



19f 2 h' 2 ^ ?// 



r 3 /i 3 Ar 2 fh 3 

f 3 7f , h ,3 



32/ 2 /i 4 r/i 5 



/2 (•// 



6r 2 fh Qrf 2 h 2 
fh" f' 2 h" 



\2f 3 h 2 rfh 3 8f 2 h 3 



9f'tif" I9ti z f 



r 2 h 3 
3h'f 



2rfh 3 
h!" 



+ 



+ 



Aph 3 
f'h 



IZh'h" + Uf'h'h" 



2rh 4 



8fh 4 



f'h" + 2f" 



8fh 4 

rr 



Afh 5 

3f' 2 
8f 2 h 2 



fh 3 rfh 2 2f 2 h 2 



III £1, 
+ ' 



2fh 3 rh 3 Afh 3 2fh 



(7.6a) 



(T e 6 



1 



240vr 2 
7/ 



7k 4 
12/ 2 



/3 



+ 



5k 2 
6r 2 / 

f 4 



5k 2 



+ 



Qr 2 fh 
f'h' 



5k 2 f 5k 2 f' 2 



+ 



2/' 



5/ 



/2 



5/ 



/2 



6r/ 2 /i 24/ 3 /i r 3 fh 2 8r 2 f 2 h 2 



fh! f' 6 h' 7f'h 



/3| 



P'A' 2 



24r/ 3 /i 2 192/ 4 /i 2 2r 2 //i 3 4r/ 2 /i 3 



+ 



2Ar 2 f 2 h 
2f" 



+ 



16/3/i 3 

r/ 2 /i 2 8/ 3 /i 2 ' rfh 3 Aph 3 8f 2 h 2 rfh 3 



2/7" ft" , W W / 



w2 



+ 



4r//t 4 
'V 



/ 



7f' 2 h' 2 



32f 2 h 4 r 2 //i 2 



(T/)o 

/' 



1 



240vr 2 
3/ 



7k 4 
12f 



5k 2 /' 5k 2 / 



2x/2 



+ 



/3 



r 3 fh 2 lQrf 3 h 2 



+ 



31fV 2 , 7/'/i 
32/ 2 /i 4 
19h' 2 f" 



+ 



+ 



8//i 4 
3h'f" 



+ 



4//i 5 

7/ //2 
8f 2 h 2 
f'h'" 



+ 



12r/ 2 /i 12/3/i 

29/ /4 3/W 

2r 2 //i 3 ~ 

7/7" 



hn 2 h' 
~ \2rfh 2 
7f' 2 h' 



+ 



8f 2 h 3 
hn 2 f'h' 



rfh 2 



Af 2 h 2 _ 
(7.6b) 
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6Aph 2 
f" 



+ 



r 2 fh 2 
f'h" 



+ 



lQrf 2 h 3 
73f' 2 f" 



+ 



5k 2 / 

24/ 2 /i 2 12/ 2 /i 

73f' 3 h' 3f'h' 2 



\2rf 2 h 2 A8f 3 h 2 
3f' 2 h" 13f'h'h" 



+ 



Arfh 3 8f 2 h 3 8fh 4 



+ 



96/3/i 3 
11 /i'f 
Arfh 3 
f'h" 



Arfh 4 
V'h'f" 
f 2 h 3 
3/"' , 3/7 



+ 



fh 3 2rfh 2 Af 2 h 2 
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2fh 3 Afh 3 2fh 2 



(7.6c) 
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[Tfi )iog 

\TfM )log 



1 



4vr 2 ! 
1 



£V In 



1/2' 



4^* 



In 



2A J 



+ C 



K 



if T = 0, 
if IV 0. 



(7.6d) 
(7.6e) 



In the final pair of expressions, // is a renormalization parameter, A is the previously men- 
tioned infrared cutoff (necessary only in the zero temperature case), C is Euler's Constant, 
and 



*V _ Jo {R PI mt u R pt - \R pT R pT g^ - -j^R [R fJ v - \Rgf 

I J_/D V\ P 1_ D V 1_ T3 P n V 

' 20^ P >',P 60 >P 120 'P Up ' 

The zero temperature case is obtained by setting k = in the above equations. 
Our results for (T A1 ' / ) num are: 



(7.7) 



r 2 f 

,2 fll 



uj 2 A 1 



r 2 u> 3 



(Tr r )i 



+ 

4 
+ 



f"\ 



Ufh 



f 

+ / 



/ 1 1 rf _ 5r 2 f 2 rh' _ r 2 f'ti 
U \U~ Uh + 6fh~ 48fh + Y2h? ~ 24fh 2 

(7.8a) 



d\i 
Z7T f 

dfi 

uJ 



^A 2 f' 2 

uj Ai + 



r 2 u 3 



1 1 



r ji r 2j 



2fi2 



12 Uh Ufh 48/ 2 /i 



Ur 



(7.8b) 



(Tq )num (-^^}num 2 J j 



00 



r 2 f'ti r 2 f" 



24fh 2 Ufh 



A 2 f' 2 2r 2 u 3 ( rf r 2 f 2 rh' 
~ + 3/ + U \Ufh ~ Ufh + Uh 2 
dfi 



(7.8c) 



where it is understood that, in the zero temperature case, there is an infrared cutoff A in 
the final integral of each equation. In these equations, the quantities A\ and A 2 are given 
in Eqs. 



VIII. DISCUSSION 

A method for computing the stress-energy tensor for massless spin | fields in static 
spherically symmetric spacetimes has been presented. The full renormalized stress-energy 
tensor can be computed either numerically or analytically in a particular spacetime geometry 
by first solving the radial mode equations (|2.37aJ -b) using the correct boundary conditions 
for the problem at hand. The result is then substituted into Eqs. ( |6.5| ) and ( |7.8aj — c), and 



the mode sums are computed. Then the analytic approximation is computed for the given 
spacetime geometry and added to the result. 

For an analytic computation this is all that is necessary. However for numerical com- 
putations there can be a problem with getting the mode sums to converge to a sufficient 
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accuracy since, from a practical point of view, one must truncate both the mode sum in £ 
and the integral (or sum) over u at large values of these quantities. Fewer modes can be 
used if one makes use of the full WKB approximation for those modes with large values of 
£ and/or u>. One way to do this would be simply to use it for modes beyond the cutoffs. 
Another is to use the full WKB approximation for the radial mode functions for all values of 
u and £ in Eqs. ( |7.8a| — c). We shall call the resulting quantity (T /li/ )wKBfin- It can be added 
and subtracted from the rest of the stress-energy tensor with the result that 

(2^L«//ren (T^y^j modes WKBfin yTfiu) analytic 

modes unrcn 

(TfjJ) WKBfin — (Tfj, u ) WKB ~ (^/W WKBdiv • (8-1) 

The quantity (T /Ji/ ) WK Bfin can sometimes be computed analytically but often must be 
computed numerically. However, even if it must be computed numerically, one can go 
beyond the necessary numerical cutoff by expanding the WKB approximation for the modes 
with very large values of £ and ui in inverse powers of these quantities. Then for each power, 
the series in £ can be summed or the integral (or series) in u can be computed analytically. 
Procedures for doing this are described in |f2T] . 

It is also possible to use (T Atjy ) ana i ytic as an analytic approximation for the stress-energy 
tensor. Examination of Eq. (|7.6a| -e) shows that our result agrees with the analytic ap- 



proximation derived by Frolov and Zel'nikov [20] when their expression is evaluated for a 



static spherically symmetric spacetime and their undetermined coefficients have the values 
= = and qf^ = j4t. Thus it also agrees with that of Brown, Ottewill, and Page |7| 
in Schwarzschild spacetime. 

Work is in progress to generalize our method of computing (T^) to the case of a massive 
spin g held in a general static spherically symmetric spacetime. We also plan to use our 
method to numerically compute the stress-energy tensor for both massless and massive spin 
\ fields in Schwarzschild and Reissner-Nordstrom spacetimes. 
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APPENDIX A: COMPARISON WITH THE USUAL FLAT SPACE GREEN'S 

FUNCTION 

In this appendix we show that in flat space, our zero temperature Green's function Se is 
equal to one derived by letting the Euclidean Dirac operator act on the usual scalar Green's 
function Q& multiplied by the unit matrix. We do this by defining a new Green's function 
G via the equation 

S E {x,xi) = - lE V*G E {x,x>) (Al) 
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and showing that 

G E (x,x') = Q E (x,x')l 
In flat space Ge can be written in the form 

/ A q > p 



Ge{x, x') = — 



with 



2irrr' 



1 m=—j 
1 = 2 



r—r ) ( XJ ' u} '^ m 
J\ o 



A q, P I + 6{r'-r) (p «-> g) 



(A2) 



(A3) 



To verify that this gives the correct Se in flat space one can first substitute Eq. (|A3| 
Eq. (|Alj) with the result that 



into 



-l E V E u G E {x,x>) 



ia (d r — -a ■ L 



(d r -}a-L 



—i 



2ixrr' 



Or 

9(r—r') x 



G E (x,x') 



j=l/2 m=-j 



uj,j,m 



-G (d r — -(rr 



1 ' - • L + 1) ) A™ jin 



a(d r -l(a-L + l))A^ m 



10 ui ,j,m 



+9(r'—r) (p <-> g) 
Next Eqs. (|2.34|) can be used to show that 



(A5) 



a [ d r — -(a ■ L + 1)) A^ m = (d r - ( -^) F^(r)F^(r')^_(6, <f>)*f + (9', 0') 

2 ' 1 Gl^Gl^r')^, 0)*™!(0', 0') . (A6) 



- & + 



r 



Finally through the use of the first order equations ( |2.37a|) and ( |2.37b| ) one obtains Eq. ( ft.2D . 
The usual form of the Euclidean Green's function in spherical coordinates in flat space 



is 



Qe{x,x') 



~ — E(2^ + l)kr-r')<,(r)< £ (r') 



+ e(r'-r)ulJr)ul e (r') P e (co S1 ) (A7) 



where cos 7 is defined in Eq. ( |3.10| ). The radial mode functions u Wj e satisfy the second-order 
differential equations 



+ 



-u ul t = . 



(A8) 
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They also satisfy the Wronskian condition 



u. 



(A9) 



Here u^ e is finite at r = and is finite at r = oo. 

To show that Ge(x,x') = Qe(x,x')1, first note that the mode functions u u> i can be 
chosen to satisfy the recursion relations 



UJU, 



d r u p 



-U 



P 



UJU, 



-U 



(A10) 



One can then easily show that both the differential equation and the Wronskian will be 
satisfied for u u t if they are satisfied for u u t-\. If we choose 



F p ■ 
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U . 1 



^Vj — M . , 1 
9 



F q ■ 



-U 



u,3- 



1 5 



(All) 



then the F's and G's will also satisfy Eqs. ( |2.37a| ) and (|2.37b| ), and the appropriate boundary 
conditions. 

Next note that substitution of Eqs. ( |A11|) and ( p.36|) into Eq. (|A3|) gives 



/c 
-I 



2nrr' 



e+7 



0(r-r')\Y: E 



-£=0 



1 L 

~2 



1 ^ m4 (M)C f (^iT 



Y; , ^{e,<p)Y i m+ ~ 2 {e',^y 
i 

2 



i 



2£+l 

DO 



( e 4^Y e m -HeA)YrHe>,<pr 



+ E E <i(r)ul/r')x 

£=1 I 1 



V - V( ^iTr m V +l (0, ^T 1-1 (c, 0')* ^K"" - < o )5-;'"" 7 ( ^. <>' ) 

+ 6{r'— r)< p ^ q 



(A12) 



The off-diagonal terms can be shown to cancel. Shifting the sums over m by ±| can then 
be shown to give 



G E (x,x') = 1® J 



oo 2irrr' 

Y e m {6,<j))Y e m {6',<j)'y 
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Using the identity flOl) one can rewrite this expression as 
Ge(x, x') = 1 

j-oo swrr — 

+ 9(r'-r)ulAr)ulJr') P e (cos^) 



a 2 ; E( 2 ^ + 1) *(r-T>Wr)«k(rO 



(A14) 



Comparison of Eqs. (|A14| ) and (|A7|) shows that Eq. (|A2|) is correct. 

Since in flat space we can obtain from Eq. it is clear that our expression for Se 

is equal to the usual flat space Green's function obtained from the relationship 



S E {x,x')=- lE V* (Ge (x,x')1 



(A15) 



APPENDIX B: BI-SPINOR OF PARALLEL TRANSPORT 

The bi-spinor of parallel transport I(x, x') is required to parallel transport a spinor at x' 
to a spinor at x P^|J?T|Jd^ ] 

ip{x)\\ — I(x,x')if)(x') . (Bl) 

From this definition it is clear that I(x,x) = 1. Furthermore, if you parallel transport ip(x') 
from x' to x, and then transport it back along the same geodesic, it should return to its 
original value, so 

I(x',x)I(x,x f ) = 1 . (B2) 

The bi-spinor can also be used to parallel transport ijj(x), namely, 

^( x )n =^( x ')7(x,x') . (B3) 

Since %fj{x)ip{x) is a scalar quantity, we expect it to be unchanged by parallel transport, so 
that 

i>{x') if)(x') = ^(x)\\ 4>{x)\\ = if)(x')I(x, x') I(x, x')if)(x') . (B4) 
Comparison with Eq. ( B2|) shows that 

7(x,x') = I(x',x) (B5) 

and therefore 

ip(x)\\ = tp(x')I(x' ,x) . (B6) 

For an object of arbitrary spin to be parallel transported along a geodesic, the inner 
product of its covariant derivative and the tangent vector of the geodesice must be zero. 
The tangent vector is equal to <r M , which is defined in Eq. (p.2|). Thus for a spin | field 

^V^(:r)|| = 0. (B7) 
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Using Eq. (pl|), one finds that J§,|§ 

a^V^z, x') = 



:bs) 



In this paper we are renormalizing by splitting the points in the time direction as shown 
explicitly in Eq. ( [5.51 ). Thus we seek an expansion for I(x, x') in powers of e = t—t'. However, 
Eq. (p8|) contains derivatives of I(x, x'), so we cannot simply assume an expansion in powers 
of e and substitute into Eq. ( |B8| ). Instead it is necessary to consider an expansion for I(x, x') 
that is valid for an arbitrary separation of the points. Such an expansion is of the form 



I(x, x') 



1 + I a (x)e a + i/ Q/3 (x)eV + I/^^eV^ + ^I a p lS {x)e^ <e> 'e 5 + ■■■ (B9) 



4! 



with e a = x° 



x 



and each coefficient I a /3... symmetric in its indices. Substitution of Eq. 



D into Eq. ([B8]) and use of Eq. ( gig ) yields 



^ [(r A 



. 



(BIO) 



One can next expand <r M in powers of e a . The procedure for doing this is described in |6 
and also in Appendix B of [ STJj , where a misprint of |64| is corrected. Q One finds 



with 



2\ b <4> 



' r \ e * e P + 



3| a/?7 



r)e°e /3 e 7 + 



4! b a/3 1 8 



V)e"e /3 e 7 e <5 



(Bll) 



s a 

S af3 
fj, 

S a/3f 



2a P 
1 a/3 

1 a/3 1 p7 "r 1 a/3,7 



,5-"- /3 7 p + Tp^R^aSp + TpgR^ 



a-yp 



_pP p/^ _ -pp -pp 

1 a/3,6 1 p7 1 a/3 1 p7,<5 



pp p/i 

1 a/3,7 1 p<5 



pP pr p£t 
1 a/3 1 P7 1 t8 



1 a/3,75 



r 75 (-R M /3pa 



ap/3j 



Substituting Eqs. (|B|) and ( |BTlD into Eq. (|BTo|) results in an equation that contains terms 
proportional to e a , e Q e^, ... . Setting the coefficients of these terms to zero results in a 
series of equations that can be solved for I a , I a/3 , and so forth. Since we ultimately want the 
points to be split in the time direction we only need the quantities I t , Itt, Ittt, and I t ttt in Eq. 
( |B9|) . However to obtain expressions for these it is also necessary to additionally compute 
I r , I tr , I rr , htr, hrr, an d I rrr . We find for these latter terms that 



7°7 





f 



J-tr 
htr 
Itrr 



f 3 f' 2 h' 



+ 



16/t 8 
f'f" 



32/ 2 /i 32fh 2 

' f 3 



Wfh 
\\f' 2 h' 
PW 2 \~l2f 2 ~ 192 fh 



7°7 



11 fh' 2 f'f" 7h'f" fh" 



192K 1 



6/ 



96h 



2Ah 



r 

12 



(B12) 



5 The equation at the top of page 109 in Ref. |34|] contains a sign error in the second term. 
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The terms that survive in Eq. ( P9|) when the points are split in the time direction are 

/'7°7 



It 
Itt 



Af l / 2 h 1 / 2 



16fh 

7°7 



( r 3 r 2 h' rr 



J «* fV2 h i/2 [ Q 4fh + 32h 2 16h 
7 f 4 f' 3 ti f' 2 f" 



itttt 



256 f 2 h 2 32fh 3 I6fh 2 
Thus the bi-spinor of parallel transport, in the time direction, is given by 



(B13) 



I(x, x') = l + I t e + \l tt e 2 + ll at e 3 + i W 4 + 0(e 5 ) . (B14) 

Similarly, 

I(x', x) = l-I t e + \l tt t 2 - \l m e 3 + l 4 I tta e 4 + 0(e 5 ) . (B15) 



APPENDIX C: DERIVATION OF AN EXPRESSION FOR THE 
UNRENORMALIZED STRESS-ENERGY TENSOR IN TERMS OF THE 
HADAMARD GREEN'S FUNCTION 



In this appendix, we use a method similar to that followed by Christensen [51], to derive 
an expression for (T^) unren in terms of the spinor Hadamard Green's function S^'(x,x'). 
However, we work with a complex field instead of the real one treated by Christensen, and 
we use different conventions for the Dirac matrices. 

The action for a massless, complex spin | field can be written |5? 



S 



d 4 x, 



' 9 2 



tfj{x) ^V^{x) - (V„V(aO) l^i x 



(CI) 



The stress-energy tensor for the classical spin | field is obtained as the variation of this 
action with respect to the metric 



5S 



- {^{x) 1{i y v) ^{x) - {V { ^{x)) lv) i){x)} . (C2) 



Using the anticommutivity of spin | fields and the cyclic property of the trace, one can show 
that 



% / r "I r 

T^(x) = -Tr (7^ \V u) il)(x),i)(x)\ + 7^ [V u )i[)(x),i/;(x 



(C3) 



Next the points can be split symmetrically, so that 
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Tp,(x) = | lim Tr {7^ (v„) ip(x),ip(x') ) I(x',x) + ^I(x,x')V l/ ) ip(x),ip(x' 

+ J(„9») X ' (Vv [^(x')^(x)]) I{x\x)+ 1{ll g u) x 'l{x,x')Vx> [^{x')^{x)]} . (C4) 

Here we have used the bi-vector of parallel transport g^' and the bi-spinors of parallel 
transport I(x, x') and I(x', x) to parallel transport the part of the expressions that are 
spinors and vectors at the point x' to the point x along the shortest geodecie connecting 
these points. 

When the fermion field ip is quantized we reinterpret ( |C4j ) as an equation for the stress- 
energy operator. One can then compute an expectation value of T^ u (x) in Eq. flC4]) with 
respect to whatever state the fields are in. It is possible to write (T^x)) in terms of the 
Hadamard Green's function which is defined as 



S^\x,x') = ( ip(x),ip(x') 



(C5) 



We find 



(T M „> = § hm Tr { 7(m (y v) S^ (x, x')) I{x' , x) - 1(ll I(x, x')V u) S^ {x\ x) 

-7(^) V (Vx>SW{x,x>)) I(x , ,x)+^ fl g u) y I(x,x')VyS^(x , ,x)} . (C6) 



To simplify Eq. 



, first note that 

S (1 \x,x') = S w (x',x) 



(C7) 



which can easily be seen from Eq. (|C5[) . Since SW{x,x') is formed from the tensor product 
of if)(x) and ip(x'), 



V^ 1 ) (x, x') = d u S^ (x, x') + r^ 1 ) (x, x') = d v S w (x', x) + S {1) (x', x)Y v . 



(C8) 



Using Eq. (|2.13| ) and Eq. (|2.5|) , together with the identity 7 M = 7^, it is easy to show that 
r,, = — r„ so that 



V v S^{x,x') = d„S {1) (x',x) - S^{x',x)T v . 



(C9) 



Since S^(x',x) is formed from the tensor product of ip(x') and "0(x), and the covariant 
derivative of a barred spinor is 

V^(x) = d v ^(x) - ${x)T v , (CIO) 
we identify the right hand of Eq. flC9|) as simply a covariant derivative, so that 

V u SW(x,x>) =V u S {1 \x',x) (Cll) 
Barring this equation, and swapping x <-> x', we see that 



V^S , ( 1 )(x,x') = V^V^) 
Substituting these two equations into Eq. (|C6l ) we see that 



(C12) 
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(T^) = § lim Tr{ 7(/ , Cv v) S^\x, x')I(x', x) - I{x, x')V v) SW(x, x>) 

- g v) x 'VxSM (x, x')I(x', x) + ^) A 7(x, ^Va^W^,^)) } (C13) 

Using the cyclic property of the trace, together with Eq. (P5|) , this can be written 

^ Tr (y {tl V u) SM(x,x')I(x',x)) - Tr ( 1( ,V u) S^(x,x')I(x' ,x)) 
Tr U„g v) x 'Vy)S w (x, x')I(x', x)) + Tr ( 7 (^) A ' Vx>)SW(x, x')I(x', xj)} (C14) 



(T^) = | lim 



One can then use the result for an arbitrary matrix A that Tr A = Tr [A]*, to write 



(T/jJ) = -\ lim Im{Tr [ 7(m (V^^rr, x') - g v) x 'VvS<»(x, x'j) I(x',x)]} (C15) 



APPENDIX D: RELATIONSHIP BETWEEN AND S E 



In this section, we establish a relationship between the spinor Hadamard Green's function 
and the Euclidean Green's function Se- This is done by first finding the relationship 
between and the Feynman Green's function Sf and its charge conjugate S F . Then a 
relationship is established between Sf and Se- 

To establish a relationship between S^ 1 ' and SV, one may start with the identity 

iS F (x, x') = x') + iS(x, x') , (Dl) 

where 

iS{x, x') = | [©(t-t') - 6>(^-*)] ({^(x), ^(O}) . (D2) 

It is possible to write S^ in terms of Sp by taking advantage of charge conjugation 
symmetry. The charge conjugation operator acts on a solution ip{x) of the Dirac equation 
by replacing it with 

^ c (x) = C$(x)) T = Cj 0T iP*(x) , (D3) 

where T means transpose and C is the charge conjugation matrix which satisfies the con- 
straints 

CjJC* = 7a and CC* = -1. (D4) 

In the Dirac representation we can choose C = i 7 2 7 °. It is then a relatively simple matter 
to show that, even in curved spacetime, if ip(x) is a solution of the Dirac equation ( p.!4|) , 
then ip c is also a solution. Furthermore, it is easy to see that (ifj c ) c = if). 

Suppose that in an arbitrary static spacetime all the positive frequency solutions are 
labeled ipk(x). Then the negative frequency solutions will be ipki x )- Hence when ip(x) is 
quantized, it can be written in terms of creation and annihilation operators as 
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ij(x)=J2[Mx)a k + r k (x)bl] ■ (D5) 
k 

Here the operators a k and b k satisfy the anticommutation relations 

{a k , a\,} = {b k , b{,} = 5 kk > , (D6) 

with all other anticommutators equal to zero. 
If Eq. (P5| ) is substitutde into Eq. (|U5|) , then 



S^(x,x') = ^2T,([Mx)a k + r k (x)bl , Mx>l>+rAx')h]) • (D7) 



k k 1 



For noninteracting fields in the vacuum state or a thermal state, the Hamiltonian and hence 
the density matrix is diagonal in the occupation number representation. Therefore, 

{a\,a k ) = hwfk , (b{h>) = <W/* , and (a k b k >) = (44') = > ( D8 ) 

where f k and f k are the occupation functions for the particle and antiparticle states respec- 
tively. Since the thermal state will, presumably, occupy the particle and antiparticle states 
equally, f k = f k . Then substituting Eq. (p6|) and Eq. (|D8|) in Eq. ( p7|) gives 



SM(x,x') = E(l"2/ fc ) {Mx)4W)l° - C 7 0T ^(x)^ T (^)7°*CV} . (D9) 

k 

The charge conjugation operator acting on a matrix operator can be defined as 

S c = C(S) T C j = -C^ T S*C ] -i° . (D10) 
Substituting Eq. (P10|) into Eq. (|D9|), and using Eq. (|D4j) gives 

S^\x,x') = E( 1 -2A){-^7 0T ^K^)^n^7°^V + 7VV fc (x)4(a;')7 7 0t 7 } • (Dll) 

k 

Then using (7°)^ = 7 and 7°7° = 1, one finds that S^ c (x, x') = S^(x,x'). In a similar 
manner it can be shown that S°(x, x') = S(x,x'). 

Letting the charge conjugation operator act on Eq. ( pi| ), and using the fact that both 
S'W and S are unchanged, gives 

- iS c F (x, x') = ±S w {x, x') - iS(x, x') . (D12) 

Combining this with Eq. (pl|) , gives 

( X; x >) = iS F (x 7 x >) - iS c F (x, x') . (D13) 

To establish a relation between Se and S F first note some results from flat space. It is well 
known that the scalar Green's functions Ge and G F are related by an analytic continuation 
HI, so it is plausible to look for a relation of the same type between S F and S F . We take 



this analytic continuation to be 

(r-r'f -> - (t-t'f + is . (D14) 
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To see that this analytic continuation is valid, consider the Euclidean Green's function Ge 
for the massless scalar field in flat space 



G E (t, x; t',x' 



47T 2 



[r - r'f + (x - x'f 



(D15) 



If it is analytically continued to the Lorentzian sector using Eq. ( pi4| ), one finds that 



Ge (it, x; it', x') 



4vr 2 



V 



(t-t'f + (x-x'f , 



—S 

47T 



- [t-t'f + [x-x'f 



(D16) 



where V denotes the principal value and use has been made of the identity from complex 
analysis that (a + ie)^ 1 = V(a^ 1 ) — m5(a), with the limit e —* + understood fl65] . From 
this equation it is easily seen that the correct relationship 



Ge (it, x; it' ,x') = iGF (t, x; t', x') 



(D17) 



is obtained |58|.p| 

One may then obtain Se from the flat space expression for Ge through the relation 

S e (t,x;t',x') = - lE Vp E (T,x;T',x') (D18) 

and Sp through the relation 

S F (t,x;t',x') =iYV ll G F (t,x;t',x') . (D19) 

Using Eqs. (|TJT7| ), ( gTg ), (|2~21| ), and ( pT23| ) it is then easy to show that in flat space 

S E (it, x- it', x') = iS F (t, x; t', x') (D20) 

where S F (it, x; it', x') is shorthand for the analytic continuation of S F (t, x; t', x') through 
the use of Eq. flDil) . 

For a general static spacetime Eq. ( |1J2(J| ) is consistent with the equations that Se and 5V 
are known to satisfy. To see this first note that the equation the Feynman Green's function 
Sf obeys is 



ij^V^Sp (t, x; t' , x" 



5 A (x — x') 



-9 



;D2i) 



while Se is a solution to Eq. (|3.1| ). If one substitutes r — it into Eq. ( |3.1| ) and uses Eqs. 
(P3) and ApD, one findsQ 



i^V^S E (it,x;it',x') 



S 4 (x — x') 



(D22) 



3 Note that in Ref. [58] the relationship r = —it rather than r = it is used. 



7 One can show, using a Wick rotation of the representation 5(t — r') = / ^e luJ ( T T \ that 
8(t-t') = -i6(t-t'). 
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Thus Se (it,x;it',x') satisfies the same equation as iSf (t,x;t' ,x') does. This provides 
strong evidence that the relationship ( P20 ) holds in a general static spacetime. 

Finally we obtain the required relation between S^ 1 ' and Se by substituting Eq. ([D2( 
into Eq. ( pi3| ) to obtain 



S {1) (t, x- t', x'\ = S E (it, x- it', x') + S C E (it, x; it', x') 



(D23) 



APPENDIX E: STRESS-ENERGY TENSOR EXPRESSED IN TERMS OF THE 

EUCLIDEAN GREEN'S FUNCTION 



In this appendix we derive explicit expressions for (T^ u ) unren and (T^) numer i c . We also 
show that both (T All/ ) nume ric and (T M y) analytic are real quantities when the analytic continuation 
( PUD is used. 

The general expression for (T MI/ ) unren in terms of Se and its derivatives is obtained by 
substituting Eq. (p23|) into Eq. ( gXp . The result is 



(T^) unrcn = -\ lim Im {Tr [ 1{fl (v„) [S E + S C E ] - g v) x ' Vy [S E + S C E }) I(x', x)] } . (El) 

To derive an explicit expression for (T A1!/ ) unren in terms of sums and integrals over the 
mode functions F u j and G u j, first substitute Eq. ( |3.11| ) into Eq. ( |ETD and compute the 
derivatives with the points split in an arbitrary direction. Next set x = x' but leave r and 
t' separated. Then the various terms in the quantity (V u )Se — g v ) X ' V \'Se} in Eq. flETp can 
be written in terms of bi-scalars multiplied by Dirac matrices. Specifically we find that 



lim Se (t, x; t', x) 



Jim 

x'^x 

Jim 

X '— *x 

lim 



Jim^ 

X 1 — >x 

Jim^ 

x '— >x 

Jim^ 
with 



S E (t, x; t , x). t \ = -7 D A + ^ 1/2fel/2 7-Pi 
S E (T,x;T',x'). t ,\ =7°D 4 + 4 , 1/a/tl/a TA 



S e (t, x: t'x') 



S e (t, x; t',x 
S e (t, x; t',x 



"7 

! 

"7 
1 

~2 



r 



D 3 + f—D 2 - 77x72-^4 + 7° ( - - — + ] Di 



r 



r 
r 



h^ 2 

'fl/2 



1 h 1 ' 2 f 



r r 4/ / 



Di 



1°D 3 + 1 



h l/2 



7°£) 1 + 7^2 



r 4/ 
(7 1 cos + 7 2 sin (fj 7 3 



(r, x: t'x'),, 



1°D 3 + 1 



h 1 / 2 



7 L) 1 + 7D 2 (7 1 cos0 + 7 2 sin0) 7 3 (E2) 



D 3 = / ^BinK^ 



D 2 = / -JJTJl^™* Ht-t')]A 2 

^j^i^J 2 cos [a;(T-r')]Ai . 



(E3) 
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The quantities A±, A 2 , and A 3 are found in Eq. ( p.5| ). In deriving these equations use has 
been made of the mode equations (|2.37a|) and ( 2.37b| ). In addition, some of the terms that 
occur in Se and its derivatives take the form 



E 



=0 L 



(e+l)G q ± F P l -iF q X G 



(E4) 



These sums are superficially divergent, but in a manner similar to that discussed in Section 
[VT| , we can add terms proportional to 8{r — r') and its derivatives to Se and its derivatives 
in Eq. (|E2j) . When the index on the first term is shifted, and the appropriate subtraction 
made, flE4|) takes the form 



E 



G q ,F P 1 -F q ,G P ! -- 



(E5) 



But by the Wronskian condition, Eq. (3/4), this vanishes. 

To derive an equation for (T^) ren one can substitute Eqs. (|E2|) into Eq. ( |E1[) and then 
compute the trace. This is more easily accomplished if use is made of the identity (7 M ) C = 
— 7 M , which holds for the Dirac representation |57]].Q The result is 



/ 1/2 Im|-(^V)p 4 -^)/i-/ r ' (i 



'■' /J + 77 ] { D i~ D i) h 



r 4/ 



tr 



+ 9 



tr 



r 

,h 1 ' 2 



(D 3 -D*)h+(g tt -g tt ' 



r 



4/1/2/ji/a { D i- D l)h 



r 



PS)* 



lm{(g"-g rr 



(D 2 -D* 2 )I 2 - g* ,l ?l(D A -Dl)l\ 

1 kl/2 1 f ) /> 1/2 Pi-£>!) / 2 - /^ 1/2 (D 4 -Dl) I 2 



4/, 



f-^) h. (D 3 -D* 3 ) h + (g rr +g rr ') ~ (D 2 -D* 2 ) h 
(g rr +g rr ') (Da-DI) h - g^-J^ (D l -D* 1 ) A J 



OA 



Im ( J D 3 - J D 3 *)/ 2 + 1 



^1/2 



{D 1 -D* l )h-{D 2 -Dl)I 1 \ . 



(E6) 



Since in each expression, we have terms of the form (D n —D^), and we are taking the 
imaginary part, it is sufficient to replace these with 2D n . This gives the expressions in 
Eq. 

In principle one may solve the mode equations (|2.37a ) and ( |2.37b|) exactly, compute the 
sums and integrals in Eqs. ([E2D, use the analytic continuation ( pi4| ), subtract off the point 
splitting counterterms displayed in Ref. [^TJ and Eq. (|7.2D , and take the limit t' — > t to 
obtain a renormalized expression for the stress-energy tensor. However, it is seldom possible 



This identity actually holds for any representation of the Dirac matrices for which = C 
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to solve the mode equations analytically and to analytically compute the sums and integrals 
in Eqs. (|E2|) . What can be done instead is to use the WKB approximation for the radial 
modes to determine the divergence structure of (T^ u ) unren . The procedure for doing this has 
been discussed in Section VII. The result is then substituted into Eq. (|E6|) to obtain the 
quantity (T M i/)wKBdiv Then 

unren unren (X/^WKBdiv) + (7ju,)wKBdiv} 

= hm {(^)nrnn + (^)wKBdiv} • (E7) 

Note that since (T^wKBdiv contains all of the divergences of (T^ u ) nnren , the quantity (T^ u ) nnm 
is explicitly finite in the limit r' — > r. As a result one may set r = r' before computing the 
mode sums and integrals in this quantity. Our explicit expressions for (T^ u ) num are given in 
Eqs. (|7T8a|-c). 



The quantity (T^wKBdiv can be computed analytically. The resulting expression con- 
tains terms of the form (r — r')~ n with n = 4,3,2,1,0 as well as terms proportional to 
In (r — r') . For even powers of n the analytic continuation takes the form 



1 ' V '— 77 I iirS 



(r-r') n -{t-t') n + ie \{t-f 
For odd powers of n one may multiply and divide by r — r': 

1 t-t' 



(t-t') n . (E8) 



r'T (r - r') n+1 



(E9) 



The numerator trivially contains no poles and therefore can be analytically continued using 
t = it, while the remainder of the expression is treated in exactly the same manner as in 
Eq. (]E8|). Upon subtracting the charge conjugate term, the terms in (T^ u ) WKBdiv that are 
proportional to delta functions cancel. For the logarithmic terms 

In (t-t') = \ In [- (t-t') 2 + ie]=\ [in | (t-t') 2 I + in] (E10) 

so that again upon subtracting the charge conjugate, the imaginary part vanishes. As a 
result, we find that (X^wKBdiv is real, and therefore that the difference (T Ail/ ) analytic = 
linv^i ((r^) wKBd iv - (T^u)vs) is both finite and real. Our expressions for (T^) ana i ytic are 
given in Eq. (|7.6a| -e). 



APPENDIX F: COMPUTATION OF THE MODE SUMS IN (T^) W kbdiv 



Here we compute the mode sums or integrals necessary to evaluate the quantity 
(^V)wKBdiv for non-zero and zero temperature states, respectively. For the thermal case, 



Eq. (1.442.2) of Reference p6| gives 



n=l/2 



cos (riKe-r 



riK 



2 



cot" 



K€ 7 
~4~ 




K 2 el 



48 



7/t 4 e 



4^4 



23040 



7> + 0(4) , (Fl) 
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where the sum is understood to be over half- integers n = |, |, . . .. For the sums and integrals 
we are interested in, e T = r — r' . Using Eq. (|F1|) and taking derivatives with respect to e T , 
we find that 

oo -y 

K £ sin (n/ce T ) = h 0(e T ) 

n=l/2 Cr 

00 Ik 2 

K E nK COS ( nKe r) = J + OA + 

n=l/2 C t ^ 4 



00 2 

k (n/t) 2 sin inner) = — - + 0(e T ) 

71=1/2 ' 

00 fi 7k 4 

k £ M 3 cos (nKe r ) = -j - — + 0(e 2 ) . (F2) 

n=l/2 e r JbU 



In the zero temperature case the mode sums are the same as in It is necessary to 

introduce an infrared cutoff A for the first integral. Then 

1,00 duj 
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/■oo 
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JO 

/■°° 2 

/ u; 2 da; sin(c<je T ) = — - 

Jo e T 

rco g 

/ cu 3 dcu cos (ue T ) = — , (F3) 
Jo el 
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where C is Euler's constant. 
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